土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.      

高熵合金的力学性能及变形行为研究进展

高熵合金是近年来提出的一种新的合金设计理念,打破了一般合金中以1种或2种元素为主,辅以极少量其他元素来改善合金性能的传统思想,由多种元素以等原子或近似等原子比混合后形成具有独特原子结构特征的单一固溶体合金.高熵合金的多主元特性使其在变形过程中表现出多重机制(包括位错机制、形变孪生、相变等)的协同,因而高熵合金已经展示了优异的力学性能,如高强、高硬、高塑性、抗高温软化、抗辐照、耐磨等,被认为是最具有应用潜力的新型高性能金属结构材料,已经成为国际固体力学和材料科学领域研究的热点.本文首先介绍了高熵合金独特的结构特征,即具有短程有序结构和严重的晶格畸变;随后对近年来针对不同类型高熵合金(包括具有面心立方相、体心立方相、密排六方相、多相以及亚稳态高熵合金)力学性能、变形行为方面的研究成果,特别是强韧化机制以及相关的原子尺度模拟,进行了较为系统的综述;最后强调了高熵合金未来研究中所面临的一些主要问题和挑战,并对其研究进行了展望.     

基于LBM方法的风力发电机尾流结构仿真

针对传统CFD数值计算方法难以实现风力机动态旋转及其旋转状态下的流固耦合计算,本文结合格子玻尔兹曼(LBM)方法易于处理动态复杂边界的特点及大涡模拟(LES)方法在非稳态涡流结构捕捉上的优势,采用LBM-LES联合方法进行三维风力发电机整机气动性能及尾流结构仿真研究,同时采用尺度自适应方法对尾涡结构进行跟踪和精细化计算。针对NREL PhaseⅥ型试验机进行模拟,得到了与实验结果吻合的流动形态及尾流结构演变规律,分析了尾流区速度演变规律并对比了不同亚格子湍流模型对计算结果的影响.     

碳离子在不同材料叶片的多叶光栅上产生的次级粒子研究

在碳离子放射治疗中，碳离子束与治疗头设备和患者身体相互作用产生的次级粒子可以到达患者体内的许多区域，在产生的次级粒子中以中子和$\gamma $射线的产额为最大。在不影响束流配送功能的情况下，减少碳离子放疗中产生的次级中子和$\gamma $射线对于降低放疗后出现的正常组织并发症及二次肿瘤风险有着非常重要的意义。本文利用蒙特卡罗(Monte Carlo)方法模拟计算了被动式束流配送系统下，400 MeV/u碳离子束照射到由不同材料叶片构成的多叶光栅(MLC)形成典型的10 cm×10 cm方形射野时，在水模体中产生的次级中子和$\gamma $射线所沉积的剂量及空间分布等。模拟结果显示:碳离子束通过MLC形成射野后在水模体中产生的次级中子主要分布在水模体的入射端，次级$\gamma $射线较为均匀的分布在整个水模体内，且较多分布在具有展宽Bragg峰(SOBP)射野在水模体中贯穿时的坪区。对于MLC叶片材料的选择，则需根据实际情况对叶片厚度以及次级粒子当量剂量的要求来确定。本文通过模拟研究不同MLC叶片材料产生次级粒子的情况，为被动式束流配送系统中MLC叶片及其他元件的材料选择提供了科学依据。     

三维格子Boltzmann方法研究上升双气泡的运动特性

应用格子Boltzmann三维模型，对双气泡在静水中的运动进行数值研究.采用八点差分和十八点差分格式分别求解一阶▽φ和二阶▽²φ可以有效避免气液密度比过大造成的数值不稳定问题.结果表明：当两个相同直径的气泡在上升时，位置靠上的气泡形状变化像单气泡上升一样，而位置靠下的气泡会受到前一个气泡尾迹的影响，并有很明显的形状变化.当两个气泡直径不同时，不管初始位置如何，大气泡总会对小气泡造成强烈的影响.     

磁化套筒惯性聚变初步探索研究与实验可行性分析北大核心CSCD

磁化套筒惯性聚变(MagLIF)构型可充分利用现有大型脉冲功率驱动装置,如聚龙一号等。基于磁流体力学方程组和1∶1比例氘氚(DT)混合燃料聚变模型,开发了零维MagLIF数值模拟程序并进行了初步探索研究。计算结果表明初始负载参数(如轴向磁场强度,预加热温度、时刻,负载半径等)与聚变产额之间有着密切的联系,在给定条件下,可依据计算给出的定性关系进行负载优化设计。值得注意的是,根据计算结果,即使在理想条件下,氘氚燃料要实现能量收支平衡,则驱动器的电流必须大于21.2 MA。这意味着聚龙一号装置(10 MA)无法开展集成化的MagLIF实验,进一步的校验计算验证了上述观点,并在此基础上提出铝套筒分解实验的建议和负载设计参数。所取得的计算结果有利于加深对MagLIF套筒压缩阶段物理过程的认知和理解。     

真空条件下LiCl溶液蒸发量随浓度变化规律研究

除湿溶液再生是除湿空调系统的核心过程,为了提高除湿溶液再生效率,以蒸发量作为评价标准,在不同的真空度下对不同浓度LiCl溶液进行再生。对不同再生条件下,溶液的再生情况进行数学模拟及实验研究,并将模拟结果与实验结果进行对比分析,结果表明,实验结果同模拟结果具有相同的变化趋势。在6 kPa和8 kPa的压强下,LiCl溶液的蒸发量随浓度的增大先上升后下降,而在1 kPa的压强下,LiCl溶液的蒸发量随浓度的增大呈现出先下降再升高接着又下降的趋势。     

行星式MOCVD反应器进口结构对AlN生长的气相反应和生长速率的影响

针对行星式MOCVD(Metal Organic Chemical Vapor Deposition)反应器进口结构对AlN生长的化学反应路径和生长速率的影响进行数值模拟研究,通过改变反应器进口形式、数量以及隔板位置发现,二重进口反应器倒置进口(即Ⅲ族在下,Ⅴ族在上)时,衬底前端的含Al粒子浓度明显升高,尤其是MMAl的浓度比传统进口反应器高两个数量级,气相反应中热解路径占主导,薄膜生长速率明显提高.在倒置进口的基础上优化隔板位置,生长速率略微降低,但薄膜均匀性明显改善.当反应器进口数量从二重变为三重和五重,反应从热解路径占主导变为热解路径和加合路径共同作用,薄膜生长速率逐渐增加,而均匀性明显改善.     

肌红蛋白力致去折叠的全原子分析

位于脊椎动物骨骼肌中的肌红蛋白在生命过程中扮演着重要角色，其折叠过程依赖于血红素的绑定.本文着重于肌红蛋白在力的作用下去折叠过程的全原子统计分析.结果表明血红素不仅具有生物学上的功能，而且决定其去折叠的动力学过程.在血红素不存在的情况下，肌红蛋白的力致去折叠过程存在一个中间态，其构象不同于化学变性剂所导致的去折叠中间态，发现新的中间态.结论与相关实验结果相符，揭示了肌红蛋白力致去折叠的一般机制.     

基于广义最大熵的具有模糊输入输出的回归模型的参数估计

模糊回归是在模糊系统中建立因变量与一组自变量之间关系的重要工具,以评估模糊自变量如何影响模糊响应变量的过程。当系统中出现小样本或者非列满秩设计矩阵时,模糊最小二乘法可能得出偏误估计。本文基于文献[13]中的多元线性回归模型,利用广义最大熵方法,针对模糊输入模糊输出数据,给出线性回归模型的参数估计和算法步骤。当输入或输出数据退化为清晰值时,该估计退化为清晰输入模糊输出或者模糊输入清晰输出的回归模型参数估计。本文结合模拟数据和实例数据,将广义最大熵方法与模糊最小二乘方法、岭估计方法进行比较研究,结果显示广义最大熵方法的可行性和有效性。